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THE INFLUENCE OF RESPONSE FEEDBACK LOOPS ON THE 


LATERAL-DIRECTIONAL DYNAMICS OF A 

VARIABLE -STABILITY TRANSPORT AIRCRAFT 

By Kenneth J. Szalai 
Flight Research Center 


SUMMARY 


Several response feedback loops are analyzed to determine their effects on the 
lateral -directional dynamics of a variable -stability transport aircraft. The response 
feedback system feeds back response variables such as sideslip angle or roll rate as 
rudder or aileron commands, or both, thus altering the various transfer functions 
which describe the dynamic characteristics of the aircraft. The range of feedback gain 
for which approximate expressions are valid in describing the effect of a particular 
loop is noted. The root -locus method is used to show how the Dutch roll, roll, and 
spiral modes are influenced as a function of feedback gain. Expressions are developed 
which directly relate feedback gains to some response parameter, such as Dutch roll 
frequency, damping ratio, or roll and spiral mode time constants. The expansion to 
multiloop systems is discussed, along with limitations of the response feedback system, 
from the standpoint of operating a variable -stability aircraft. 


INTRODUCTION 


In conjunction with the design and development of a general purpose airborne simu- 
lator-*- (GPAS), studies were made at the NASA Flight Research Center of response 
feedback control loops and their influence on the GPAS, a variable -stability Lockheed 
JetStar aircraft. Although much information has been published on variable -stability 
aircraft, no single document contains a development of the fundamental concepts of 
response feedback systems along with a detailed discussion of loops that are utilized in 
variable-stability aircraft. Reference 1 is concerned primarily with showing the fea- 
sibility of installing a response feedback system (RFS) in a T-33 aircraft. The theo- 
retical treatment is limited to showing that the aerodynamic coefficients of the equations 
of motion of the T-33 can be altered to match the coefficients of an aircraft to be simu- 
lated. The report does not analyze specific loops in detail, nor does it predict how 
individual loops influence the dynamic modes of the aircraft. It is not an adequate 
source of basic information on analysis techniques. Reference 2 is typical of many 
reports which emphasize research results of using an RFS aircraft but do not present 
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much background material on the theoretical basis of the system. This particular 
reference does contain several approximations derived from the equations of motion 
which are useful in describing the effects of certain RFS loops. 

A qualitative description of a response feedback system is given in reference 3 , 
which discusses airborne simulation in general and its applications in research and 
training. Reference 4 is devoted entirely to one loop configuration and its effect on the 
spiral mode. It does contain some expressions that are useful in the description of the 
effect of several loops on the spiral mode. Reference 5 correlates root -locus analysis 
with the equivalent stability -derivative approach in the description of single-sensor 
control loops. The discussion is limited to two specific feedback loops and is not pri- 
marily concerned with variable-stability aircraft applications. 

This paper summarizes the results of an analysis of the influence of response 
feedback loops on the lateral -directional dynamics of a variable -stability transport 
aircraft. It is intended to fill a significant gap in the available literature by presenting 
fundamental information on the response feedback system, demonstrating some of the 
more successful analysis techniques, and discussing in detail the feedback loops 
commonly used in variable -stability aircraft. Somewhat unique is the range of feed- 
back gain for the JetStar for which certain classical approximations are valid in 
describing the effect of a particular loop. 


SYMBOLS 


The units used for physical quantities in this paper are given, where applicable, 
in both the International System of Units (SI) and U. S. Customary Units. Factors 
relating the two systems are included in reference 6. 

The body-axis system and sign conventions used in this paper are shown in fig- 
ure 1. 

A,B,C,D,E coefficients of fourth-order algebraic equation 
b wing span, m (ft) 

Cj) drag coefficient 

q rolling -moment coef- 

ficient 

C n yawing-moment coef- 

ficient 

Cy side -force coefficient 

aq 

VW 



Figure 1.— Body-axis system and sign conventions. 
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e 

g 

*XX» I YY* *ZZ 
X XZ 


natural log base 

acceleration of gravity, m/sec 2 (ft/sec 2 ) 

moments of inertia referred to X, Y, and Z body axes, 
respectively, kg-m 2 (slug -ft 2 ) 

product of inertia referred to X and Z body axes, kg-m 2 (slug -ft 2 ) 


j complex variable operator 

jo, designation of imaginary axis of s -plane, rad/sec (deg/sec) 

K a constant 

t = S(Sb _b_ 

Ixx2V% 


t = b _ C7 

r I xx 2V l r 


t = a^k 
p Ixx 


Cz 


ft 


H 


= a§b b_ 

Ixx 2V % 



qSb 

x xx 





L 




artificially created stability derivative; rolling acceleration due to 
bank angle 


m 


mass, kg (slugs) 


jj _ flSb _b_ _ 

N P I ZZ 2VS> 

N 0 ‘ % 

N /j = ifj h c *e 


H »a‘l 


qSb 


ZZ 




a 


N «r = ^ S>. 


N 

P 

q 


<p 


r 

S 

s 

u 

V 

V 

w 


artificially created stability derivative; yawing acceleration due to 
bank angle 

rolling angular velocity, rad /sec (deg/sec) 
pitching angular velocity, rad/sec (deg/sec) 

dynamic pressure, |pV^, N/m^ (lb/ft^) 

yawing angular velocity, rad/sec (deg/sec) 
wing area, m^ (ft^) 

Laplace transform variable 

velocity component along X body axis, m/sec (ft/sec) 
true airspeed, m/sec (ft/sec) 

velocity component along Y body axis, m/sec (ft/sec) 
velocity component along Z body axis , m/ sec (ft/ sec) 


v = qs b 

r m v 2V C Y r 
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X, Y, Z 


body-axis coordinates 
angle of attack, rad (deg) 


a 

P 


' a 



£<p 

p 


angle of sideslip, positive for velocity vector toward right, rad 
(deg) 


total aileron deflection, 6 a ^ e ^ - 6a r ight* P os ^^ ve f° r left 
aileron deflected trailing edge down, rad (deg) 


pilot aileron command, positive when commanding positive S a , rad 
(deg) 


rudder deflection, positive for trailing edge left, rad (deg) 


command to rudder actuator or aircraft transfer function, rad (deg) 


pilot rudder command, positive when commanding positive 5 r , rad 
(deg) 

damping ratio of second-order dynamic response 

-^-(s) transfer -function numerator parameter 
°a 

Dutch roll damping ratio 


9 pitch-attitude angle, rad (deg) 

p mass density of air, kg/m 3 (slugs/ft 3 ) 

<j designation of real axis of s -plane 

T r roll mode time constant, sec 

t s spiral mode time constant, sec 


<P 


roll angle, rad (deg) 


<?A 

A 

0 

w 




actuator phase lag, rad (deg) 
arbitrary phase angle, rad (deg) 
heading angle, rad (deg) 
frequency, rad/sec 

undamped natural frequency of a second -order dynamic response, 
rad/sec 

k~(s) transfer -function numerator parameter 
°a 

Dutch roll undamped natural frequency, rad/sec 


Feedback loop and transfer-function nomenclature 


1 . 


-S. feedback loop which commands aileron proportional to roll rate, and 

P Fi 

the actual value of feedback gain (e. g. , — = 5 rad/rad/sec), 
rad/rad/sec (deg/deg/sec) p 

^a feedback loop which commands aileron proportional to yaw rate, and 

r a gain value, rad/rad/sec (deg/deg/sec) 


feedback loop which commands aileron proportional to angle of side- 
slip, and a gain value , rad/rad (deg/deg) 


feedback loop which commands aileron proportional to roll angle, 
and a gain value, rad/rad (deg/deg) 

feedback loop which commands rudder proportional to yaw rate , and 
a gain value, rad/rad/sec (deg/deg/sec) 

feedback loop which commands rudder proportional to angle of side- 
slip, and a gain value, rad/rad (deg/deg) 


6r 


feedback loop which commands rudder proportional to time derivative 
of sideslip angle, and a gain value, rad /rad/sec (deg/deg/sec) 


^Values of feedback loop gains and transfer function ratios in this paper are the 
same in both the SI and U. S. Customary Units. 
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actuator transfer function 

roll rate to aileron transfer function 

magnitude of roll angle to sideslip-angle transfer function evaluated 
at the Dutch roll root 

handling -qualities parameter 


j* 

■z — (s) yaw rate to rudder transfer function 

°r 


Subscript: 

T 


trim value 


A dot over a quantity indicates the differentiation with respect to time; a double 
dot indicates the second differentiation with respect to time. 


Special notation: 
basic 

closed loop 
effective 
open loop 


relates quantity to unaugmented aircraft, no feedback gains 

with at least one feedback gain ^0 

including the augmentation of response feedback loops 

all feedback gains = 0 , which results in a basic configuration 


GENERAL DISCUSSION 


In many cases the response feedback loops effectively augment the aerodynamic 
stability derivatives of the aircraft, hence the influence of these loops can be described 
by considering the usual effect of those particular stability derivatives on aircraft 
dynamics. However, some loops create aerodynamic stability derivatives that are 
nonexistent for the basic aircraft. The effect of these loops can be interpreted by 
noting the effect of the created stability derivatives. Where the classical approxi- 
mations are not valid, root -locus diagrams and analog -computer results are used to 
describe the action of the loops on the airplane responses. These other methods are 
necessary when feedback gains are high enough to augment some aerodynamic stability 
derivatives which, because of their usual small contribution to a certain dynamic mode, 
are not considered in the usual approximations. 
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Since the response feedback system can often be described in terms of equivalent 
stability derivatives, the effects of these loops vary somewhat, depending on the char- 
acteristics of the aircraft being considered. The quantitative results in this paper were 
obtained for the GPAS from analytical studies of a Lockheed JetStar, an executive -type 
jet transport with aft -fuselage mounted engines. References to loop gain values should 
be indicative of values for other types of high-performance jet transports. Most large 
transports exhibit similar lateral -directional dynamics, which indicates similar 
transfer-function characteristics, and small variations in the basic aircraft transfer 
functions generally do not alter the dominant influence of a particular feedback loop. 
Thus, a particular loop described for the JetStar will have much the same effect on 
other aircraft in its class with some variation in the amount of feedback gain required 
to achieve certain closed-loop dynamics. 

Since the individual feedback loops utilize the conventional aircraft control sur- 
faces, the effectiveness of the loop is also a function of dynamic pressure and moments 
of inertia. Specific gains mentioned are valid only for a particular flight condition. 
Numerous references are made to flight condition, thus a code is used for convenience. 
Three parameters are written together in the following order: Mach number, weight, 
altitude (in thousands of meters (thousands of feet)). Thus, 0. 55H6. 1 denotes a Mach 
number of 0.55, heavy (17,236 kg (38,000 lb)), and 6100 meters (20,000 feet) altitude, 
and 0.23L0 denotes a Mach number of 0.23, light (10,886 kg (24,000 lb)), at sea level. 

The particular feedback loops discussed herein were chosen because they are the 
most effective for varying aircraft dynamics and are used on other variable -stability 
aircraft. Most root -locus diagrams presented are for the 0.23L0 flight condition. The 
analysis was originally made for this case because it represented a low -dynamic- 
pressure condition, and hence gave information on minimum loop effectiveness. There 
is no other significance to the selection of this flight condition. Root -locus diagrams 
for this case are representative of those throughout the flight envelope. Where signifi- 
cant differences do occur, another root -locus diagram or special comment is included. 

The response feedback loops analyzed are listed in the following table along with 
their estimated maximum gains, primary influence, and reasons for reaching a limit. 
These maximums were anticipated by the Cornell Aeronautical Laboratory on the basis 
of their experience with variable -stability aircraft and were used to predict the simu- 
lation capability of the GPAS. Actual limits can be determined only by flight test. 


SUMMARY OF RESPONSE FEEDBACK LOOPS 


Feedback 

loop 

Primary 

derivative 

changed 

Estimated maximum 
range of gain 

Primary influence 

Reasons for reaching estimated 
maximum 

6 r 

0 

N 0 

±10 rad/ rad 

Dutch roll frequency, directional stability 

Structural instability and noise 

^vane 

N /j 

±1 rad/rad/sec 

Dutch roll damping 

Noise 

8r 

r 

N r 

*4 rad/rad/sec 

Dutch roll damping 

Sufficient for the task 

<P 

L * 

±5 rad/ rad 

Roll-attitude stabilization 

Sufficient for the task 

/a 

P 

L P 

±2 rad/ rad/sec 

Roll mode time constant 

Structural instability and noise 

£a 

r 

L r 

±4 rad/rad/sec 

Spiral mode time constant 

Sufficient for the task 

6a 

0 

h 0 

±10 rad/ rad 

Roll to sideslip ratio, dihedral effect 

Noise 
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ANALYSIS OF SINGLE CONTROL LOOPS 


In the following sections, two assumptions are made concerning the feedback loops 
that are analyzed: 

1. The response sensors have no dynamics of their own. They yield a 
signal proportional to the response with no time lag. 

2. The actuators are also considered to be ideal. Actuator dynamics are 
neglected. 

These assumptions are made in order to simplify the analysis to the point where basic 
loop effects are not clouded by hardware -imposed restrictions. 

6 r 

Each feedback loop can be represented as shown in figure 2, in which the — 

loop is used as an example. The commanded rudder 6 r is a combination of the pilot 

input 5^ and the rudder commanded by 
r P /6 r \ 

the feedback loop rl — J as follows: 

6 r'% +r (v) (1 > 

6 r 

Note that the feedback gain — can be 

either positive or negative. The rudder 
deflection, then, is a function of the air- 
craft’s yaw-rate response. The aircraft 
transfer function is a mathematical model 
used to describe the input -output 
relationship between an aircraft control- 

r 

ler and an aircraft response. The *”-(s) 

°r 

transfer function is calculated from the 
equations of motion (appendix A). Sub- 
stitution of equation (1) into the equations 
of motion will eliminate 6 r and add three new terms to the equations: 


— (s) transfer function 
or p 

Aircraft" 1 
transfer 
function 


+-o 


-(s) 


■ft) 




Figure 2.— Feedback-loop configuration. 



Yg^r, augmenting 
Lg^r, augmenting 
Ng^r, augmenting 


Y r r 


L r r 


N r r 
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The *~-(s) transfer function can be calculated from these modified equations and 
6r P 6 r 

• will be a function of the feedback gain — . The feedback loop will not affect all 
three terms equally but will depend on the relative magnitudes of Yg r> Lg^, 

• and Ng . If, for example, Ng r » Lg and Ng^ » Yg r> the -p- feedback loop 
would act almost entirely as an N r augmentation loop. The equations will 


also be in terms of 5r n instead of 6 t 


It is apparent that the 


loop does affect the equations of motion and the 


; r p 


-(s) 


transfer function will also be affected. It might be asked, then, which of the other 
transfer functions are changed, and, since the transfer functions are usually presented 
in the form of one polynomial in s divided by another polynomial in s, whether the 
numerator or denominator or both will be a function of the feedback gain. The de- 
nominators of all lateral -directional transfer functions are identical. This polynomial, 
when set equal to zero, is referred to as the characteristic equation. The character- 
istic equation is altered by any feedback loop as shown in references 7 and 8. The 
effect of feedback loops on transfer -function numerators is not obvious (see, for 
example, ref. 9), but the following facts are noteworthy and can be easily proved by 
deriving particular transfer functions which fall into the categories mentioned: 


1. The closed -loop transfer function resulting from feeding back a response 
variable to one control input has the same numerator as the open-loop transfer 

r 6 r 

function. In other words , the numerator of the « — (s) transfer function with a — 

°r p r 

r 

feedback loop is identical to the numerator of the * — (s) transfer function with no 

°r p 

feedbacks, which is, of course, the open-loop transfer function. 


2, Loops which feed back some response variable to a particular control in- 
put will not alter the numerators of transfer functions relating any other response 

6 r 


variable to that same control input. 


the ^~(s) or ^-(s) transfer function numerators. 


That is, the — feedback loop will not alter 


^ 6 r 

3. It is possible to alter the numerator of a particular transfer function by 


feeding back a different response variable to a different control input. Thus, the 

Q ^r 

numerator of the *£-(s) transfer function will be altered by the — feedback loop. 

°a r B r 

It will not be altered by the feedback loop or the * feedback loop. 

On Or 


This paper draws heavily from the interpretation of dynamic stability from s -plane 
considerations. These relationships are discussed in references 7 and 8 and in ap- 
pendix B, which summarizes the application of root -locus results to dynamic stability 
analysis. 
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The 


Feedback Loop 


The — feedback loop, or yaw damper, is common in many aircraft. One con- 
venient method of analyzing this loop is to consider the two-degree-of-freedom 
approximation to the Dutch roll. By including the feedback equation in the equations of 
motion in appendix A, the characteristic equation can be written approximately, as 
derived in appendix C, as 

S 2 - K + Vp * (t) N «rl s + N /3 + N r Y 0 + (t) ( N «r Y 0 ' N P Y S r ) = 0 < 2 > 

Comparing equation (2) with the standard form of the second-order approximation to the 
Dutch roll, s^ + 2 f^co^s + u> ^ = 0, results in 


2 S i - - |n r + Yp + N5J 
» 2 = Nfl + + (^) (N 5r Y|3 - MpY^) 


It is apparent from equations ( 3) and (4) that the 


loop will affect both the 


natural frequency and the damping ratio of the Dutch roll with varying degrees of 
effectiveness. For any given aircraft, numbers can be substituted into equations (3) 

6 r 

and (4) to find how effective — is in changing fa and cua. For relatively small 
5 r» /6 A 

values of — , I— INfi,, is of the same magnitude as N v + Y« for the JetStar at mos 


values of — , jN 5 r is of the same magnitude as N r + for the JetStar at most 
flight conditions. Hence, for gain values of between 1 and 2. 5 rad/rad/sec 
(deg/deg/sec), the 2f^u>^ term is greatly changed. For the same range of gain, 
however, the (^J(^8 r Yp - term is small compared to N^+NyY^. There- 

fore, little c Ojp variation is to be expected for small ^7 gain values; thus, in 
equation (3) f^ will be primarily affected. By neglecting the 00 ^ variation and 

noting that No »N r Y^ for the JetStar as well as most transport aircraft, equation (2) 
becomes 


N r + Y p + 


Us + = 0 


Thus, the — loop is almost entirely a Dutch roll damping loop for relatively small 

gain values, on the basis of the classical two-degree-of-freedom Dutch roll approxi- 
mation. 
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The root -locus technique can be used to relate the dynamic characteristics of the 

6 r 

aircraft directly to the value of the feedback gain — . Figure 3 is typical of most Jet- 
Star flight conditions. The locus shown is for a positive yaw rate fed back as a 



6 r 

positive rudder input. The pole locations for — =0 represent the open-loop aircraft. 

6 r ^ 

For — < 2. 5 rad/rad/sec, the Dutch roll damping ratio is affected, with little change 
in ufy, as the previous discussion predicted. The Dutch roll approximation does not 


adequately describe the effect of the feedback loop for gains greater than 

5 

2.5 rad/rad/sec. Equation (5), then, is valid only for —■ < 2.5 rad/rad/sec. Above 

this value of gain the Dutch roll damping ratio is still increased, but a marked increase 
in frequency is also evident. Hardware limitations (table, page 9) are expected to 
restrict the maximum feedback gain to about ±4 rad/rad/sec for the JetStar, so for 
moderate gain values the feedback loop can be considered as a damping loop for the 
Dutch roll. The dashed line in figure 3 shows how the root -locus diagram would look if 
the feedback loop had a pure effect on the Dutch roll damping ratio. 


The — feedback loop has some other noteworthy effects on the aircraft dynamics. 

As the Dutch roll becomes more stable, the spiral mode becomes more stable. (For 
a few flight conditions, 0.75L20, for instance, the zero locations are such that spiral 

divergence results. ) The roll and spiral modes merge for — = 3. 7 rad/rad/sec, 

and for greater gain values a lateral phugoid mode exists. It should be noted that 
this mode may not be attainable in practice, since a gain of 3. 7 rad/rad/sec is 
near the expected maximum. 
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JCO 


The root locus for the opposite sign of the feedback gain is shown in figure 4. 

Since root -locus diagrams are symmetrical with respect to the real axis, only the top 

half is shown in this and succeeding dia- 
grams. For very low feedback gain levels 
the spiral mode becomes divergent. As 
the gain level is increased, the Dutch roll 
becomes divergent as a result of the de- 
crease in damping ratio to zero. The roll 
mode is virtually unaffected by the feed- 
back loop in this case. If this loop is to 
be used to decrease Dutch roll damping, 
some means must be provided to maintain 
spiral stability. 



jw 


The root locus for a higher dynamic 
pressure flight condition, 0.55H6. 1 
(0. 55H20), is presented in figure 5. The 
dynamic pressure at this flight condition 
is 9911 N/m 2 (206 lb/ft 2 ), compared with 
3735 N/m 2 (78 lb/ft 2 ) at 0. 23L0. The 
roll mode time constant is larger, hence 
the pole is nearer the imaginary axis than 
before. The Dutch roll locus is nearly a 
perfect quadrant of a circle, which means 

6 r 

that for an increased — gain an increase 

in can be realized with virtually no 

change in Dutch roll natural frequency 
for £; f) between 0 and 1. The maximum 

value of gain for which no change in 

occurs is 2. 1 rad/rad/sec, approximately 
the value for which ctty started to be 

affected at the 0.23L0 flight condition, 
even though the root locus is slightly dif- 
ferent. The first breakaway into a lateral 
phugoid occurs with an attainable feed- 
back gain, 2. 3 rad/rad/sec. However, 
the high damping ratio (£ n >0.9) makes it 

difficult to detect. The second breakaway 
occurs for a gain of 4 rad/rad/sec, which is the expected maximum gain. 



6 r 

One of the practical problems in using the — loop occurs during a steady turn. 
For the stabilizing sense of the gain, the — loop commands a rudder deflection to 

oppose any yaw rate. In a turn this means that an increased amount of rudder-pedal 
input by the pilot is required to maintain some steady yaw rate. The usual method of 

6 r 

eliminating this problem is by inserting a "washout" circuit in the — feedback loop. 
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which reduces the effectiveness of the loop at low frequencies. Hence, during a steady 
turn, the yaw-rate feedback gain is zero. 


6 r 

The — Feedback Loop 

0 


For certain wings -level maneuvers, /3 can be assumed to be approximately equal 

5 r 

to -r. It would seem, then, that the effect of the ±— loop is similar in its effect 

6 r 0 

to that of the + — loop. Comparison of the root -locus diagrams in figures 3 and 6 

5 r 

shows the similarity. An expression can be derived relating the feedback -r- to a 

Dutch roll characteristic and to aerodynamic stability derivatives of the basic air- 
craft. This relationship is shown in the following equation: 


5 r + N r + Yp 

7 ~ N S r + Y 6 r ( N r + 2 


The derivation of this equation and the conditions of validity (appendix D) illustrate a 
useful approach which results in similar expressions. Equation (6) allows the gain to 
be calculated that will yield the value of 2 £^c<j^ desired. 


The root -locus plot of figure 6 is for the feedback configuration commanding a 

6 r 

negative 6 r for a positive $ . The figure shows that the — loop alters the Dutch roll 
damping ratio with little effect on hence in equation (6) the selection of some 


2£^,o>^ implies the choice of a particular 

jw 



(Gain values in rad/rad/sec.) 


with remaining nearly constant at its 

open-loop value. Equation (6) assumes 
that ojjp does not vary, and figure 6 

indicates that this assumption is good for 
5 T* I 

— < 2. 5 rad/rad/sec. Over a fairly 
arge range of gain where is varying 
rapidly ^r-j between 0 and 2 rad/rad/se<^, 
the roll and spiral modes are not influ - 

S r 

enced as much as when — is used. 

r 

If /3 is obtained by differentiating 
/3 as measured by a vane, the loop may 
be noisy, thus limiting the feedback gain 
to relatively small values. Gains of 
much more than 1 rad/ rad/ sec are not to 

be expected with a /3 signal derived by 
differentiation; hence, equation (6) is ap- 
plicable for all practical gain values. To 
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circumvent the noisy (3 vane signal, (3 may be synthesized by measuring bank angle, 
yaw rate, and lateral acceleration and solving a simplified side-force equation (ref. 10). 

The resulting (3 will be less susceptible to gust disturbances and will allow a higher 
feedback gain to be used. 


Reversing the sign of the feedback gain results in the root -locus diagram in 
figure 7. The^-oll and spiral modes are still unaffected. The Dutch roll becomes 

divergent for -£ = 0.3 rad/rad/sec. The Dutch roll damping ratio is again altered 
[3 o r 5r 

with little change in co^. Comparison of the root-locus diagrams of the — and — 

{3 

loops (figs. 2,3, 5, and 6) shows the latter to be more nearly a pure Dutch roll 
damping loop, in that the roll and spiral modes are not influenced as significantly as 

6 r 

thev are with the — feedback 


Dutch roll 



they are with the 
loop. 


The — loop adds the fol- 

(3 

lowing three terms to the equa- 
tions of motion: 
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Figure 7.— Root locus for the ~r loop at 0.23L0. 

£ 

(Gain values in rad/rad/sec.) 


fK J 


These terms may be thought of as augmenting the three aerodynamic stability deriva- 
tives Yp, ~L'p, and Np. These stability derivatives do not appear in the equations of 

motion for the JetStar (appendix A) because of their extremely small contribution. 

6 r 

The — loop, then, can be regarded as actually creating the three derivatives, 

although they do exist for the basic JetStar. In perturbation -type maneuvers, these 
three derivatives have effects on lateral -directional dynamics similar to those of Y r , 
L r , and N r . 


Feedback Loop 


Angle of sideslip is detected and fed back as a rudder command to form the — 

loop. The two -degree -of-freedom approximation to the Dutch roll is derived (appen- 
dix C) by writing an auxiliary equation similar to equation (1) and then solving for the 
characteristic equation as follows: 
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s - N r + Yjg + 1 


s + N^ + l 


n 6 -o 


6 r /6r\ 

For small values of — , the (— 1Y term is small compared to N r + Y^. The term 
^-^N 5 r , on the other hand, is comparable to N^g, hence the Dutch roll frequency will 


be primarily affected. When 


is of the proper sign to increase au, the total system 


damping 2 will increase. The two -degree -of ^ -freedom approximation (eq. (7)) is 

generally good for flight conditions of the JetStar where trim angle of attack is small. 

If this is true, the following expression can be written relating the Dutch roll frequency 

5 r 

and the feedback gain — : 


The Np approximation for is not satisfactory for large trim angles of attack. It 

is necessary to use a more complete expression for c which includes the term 
avpLg, since this term is significant for large av^. v 
6 r 

Note that the — loop acts as an effective N^, also augmenting Y ^ and to a 
lesser degree, through the terms and g r j, respectively. 

A typical root -locus diagram 
» -i j 4 for the JetStar is shown in fig- 

\ ure 8. The locus represents the 

\ configuration in which positive 

V /3 commands negative 6 r . The 

p ~ 3 - 6 \ i 3 roll and spiral modes are altered 

very little. There is some slight 

= _ 1<5 | spiral destabilization, but even 

P ’ for very high gain the closed-loop 

| " J 2 root can be no farther into the 

1 right half of the complex plane 

Dutch rollX than the open-loop zero. For the 

JetStar, this means that the time 
"" J 1 to double amplitude for the di- 

vergent term is no less than 
40 seconds for any flight condi- 
Roll Spiral 6 r 

1q qJl i i cr tion when the -— loop is being 


= -1.5 


Dutch roll 


Spiral 


P ~ 2 '* 

for — =0 
T s 
-jl 


Figure 8.— Root locus for the loop at 0.23L0. 

(Cain values in rad/ rad.) 


tion when the - 
used alone. 


The effect on the Dutch roll is 
of primary importance, and the 
poles move in the manner predicted 
analytically, following a line of 
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nearly constant £w,um,. Since the Dutch roll mode poles are usually lightly damped for 

Y ^ /6 \ 
this class of aircraft, the small increase in caused by \JfJ Y 6 usually results in 

the locus following a line of constant The Dutch roll undamped natural frequency 

6 r 


co. 




is shown as a function of — in figure 9. 

r 


rad/sec 



/ 3 

va 


he plot is fairly linear for 

6 r I 

< 6 rad/ rad. Above this gain 

6 r 

ue, — becomes less effective in 

r 

increasing co^ in the higher gain 

Practical considerations re- 
to a maximum value not 


Figure 9.— Variation of ( O ,i. with the ~~T~ feedback 


gain at 0.23L0. 


P 



ranges. 

|6r | 

strict Hj-| 

much greater than 10 rad/rad. The 
plot indicates that the effort required 
to increase the upper limit of the gain 
from that in table I might not be 
justified in light of the small increase 
realized in ci>^. 

Reversing the sign of the feedback 
gain results in the root locus of fig- 
ure 10. The Dutch roll frequency is 
reduced, and at the same time the 
damping ratio decreases. For 
6 r I / 

>2.2 rad/ rad the Dutch roll is 


unstable. The spiral mode is sta- 
bilized. The roll and spiral modes 
merge and form an oscillatory mode. 
This mode may not be possible to 
attain in practice because of Dutch roll 
instability. 

It has been noted from figure 8 

6 r 

that the — loop would change the 

total damping of the system very little. 
Figure 10 indicates that 2£^co^ de- 
creases rapidly and goes to zero for a 
fairly low value of feedback gain. This 
indicates that the two -degree -of - 
freedom approximation is no longer 
a valid representation of the Dutch roll. From equation (7) it is noted that 


(Gain values in rad/rad.) 


N r - Y p 


-(f) 


y 6i 


( 9 ) 
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By substituting actual values for the stability derivatives , it is found that for 
6 r / 

— = 2.2 rad/rad, = 0. 16. The root locus, however, shows neutral stability. 

I— ratio is high. Reference 5 discusses conditions of 


this gain area the effective 


N 


J8 


In 


validity for the 2f^a>^ approximation. For large 
by equation (9) is not valid. 


No 


ratios the approximation given 


The primary influence of — on the aircraft dynamics, that of changing alone, 

prevails for the stabilizing sense of the gain. Since the response feedback system 
works through the aerodynamic control derivatives, the effectiveness of the loops is a 
function of dynamic pressure. This dependency is shown by writing the expression for 


uty 2 as 


“» 2 = fi[ c v(^) Cn s, 


( 10 ) 


0 r 9 

It is evident from this equation that, for a given — gain, is directly propor- 

tional to dynamic pressure for regions of the flight envelope where C n ^ and 

are constant. The ratio between the closed loop o>^ 2 and the open loop co^ 2 , 

6 r 

however, will be constant in this region for a given — gain value, as shown by 
forming the ratio " 


^ closed loop 
2 

^ open loop 




' n 6, 




( 11 ) 


The 


Feedback Loop 


This feedback loop senses roll rate and commands aileron deflection propor- 
tionately. By writing an auxiliary equation similar to equation (1) and substituting it 
into the equations of motion (appendix A), it is noted that two new terms appear: 


(^r) L 5 a P, augmenting L p p 
(~) N 6 a p > augmenting N p p 


For the transport class of aircraft, L$ is usually very much greater than so 

6a 

that the 


— loop inf luences the effective L p almost entirely. The aerodynamic sta- 


bility derivative L p is strongly related to the roll mode time constant. In fact, the 
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usual approximation, which is good for all JetStar flight conditions, is 


p (effective) 

Oa 

Hence , the — loop exercises strong and direct control over the roll mode time 
P 6^ 

constant. The — gain value must be relatively large before N_ augmentation is 
P P 

significant, compared with the strong influence on the effective L p through Lg a# 


Dutch roll 




Spiral 


The root -locus diagram for the 

6 a 

loop is shown in figure 11. The 

spiral mode remains unchanged, and 
even high gains do not noticeably modify 
Dutch roll dynamics. The roll mode is 
stabilized. By using equation (12), the 
6 a 

— gain required for a particular roll 

mode time constant can be calculated 
to be 

■— + L. 


Figure 11.— Root locus for the loop at 0.23L0. 

(Cain values in rad/rad/sec.) 


Exact solution- 


-Equation (13) solution 
(Lp variation only) 


The accuracy of equation (13) is shown 
in figure 12. The solid line represents 
the roll mode time constants resulting 

6 a 

from a variation of — obtained from 

P 

an exact solution to the equations of 
motion for the 0.55H6. 1 (0. 55H20) 
flight condition. The dashed line 
represents the roll mode time constants 


and associated 


gains obtained from 


equation (13) which assumes only Lp 

augmentation. The close agreement 
shows the almost sole dependency of 
r r on Lp and indicates that equa- 
tion (13) is a good approximation for 
the JetStar. 


— gain, rad/rad/sec 

r 

Figure 12.— Comparison of exact solution with equa- 
tion (13) solution. 


For a positive 


gain, the root 


locus of figure 13 results. The roll 
mode is now destabilized. The roll and 
spiral modes form a lateral phugoid as 


20 



jw 


1 j 2 


Dutch roll 



in 


6 a 

— increases and eventually become un- 
P 

stable. The Dutch roll frequency is 
affected slightly, but the damping ratio 
is decreased noticeably. 



-Hi 


Figure 13.— Root locus for the -p- loop at 0.23L0. 
(Gain values in rad/rad/ see.) 


For roll-mode analysis of the Jet- 
Star, the complex pole and zero pairs 
can be deleted without invalidating the 
locus on the real axis. This is equiva- 


lent to writing the open loop 
transfer function 


_e_ 


(s) 


6> = 


Ks(s 2 + 2£yUyS + a 

( s+ ^)( s+ ^)P +2 w + V) 


(14) 


and noting that for the JetStar the two second-order terms are close to cancellation. 


This fact can be interpreted in terms of the handling-qualities ratio — r . For the 


co 


>P 

JetStar and similarly configured transport aircraft, this ratio is approximately 1. The 
complex pole and zero will be very close to each other, depending on how close the 


ratio 


•ip 


is to 1. Thus, the Dutch roll poles will usually close on the complex zeros, 


not influencing the locus on the real axis. 


The 


^a 

£ 


Feedback Loop 


This feedback loop commands aileron deflection proportional to angle of sideslip. 

will be modified through the 




For the JetStar and similar aircraft, Lg a is usually much greater in magnitude than 

Ng a , hence, augmentation will dominate. The influence of this loop varies more 

from one flight condition to another than the loops previously discussed but generally 
has a strong influence on the spiral mode and a noticeable effect on the Dutch roll mode. 
The action of this loop in providing a rolling moment due to sideslip is recognized as 
the dihedral effect. When Lg is increased negatively, the spiral mode is generally 

stabilized and the Dutch roll mode destabilized. 


The aerodynamic stability derivatives and N^> 
following new terms: 


{t) l ^ 


augmenting 




Ng J3, augmenting 
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loop is shown in figure 14. The 


jw 


Oo 

A representative root -locus diagram for the - — 

exact shape of the locus for different flight conditions varies somewhat because of the 

drastic variation in S~(s) transfer 

°a 

function zeros. The figure shows the 
spiral mode being stabilized and the 
Dutch roll being destabilized, as is pre- 
dicted from li£ considerations alone. 

The Dutch roll mode usually goes un- 

6 a 

stable for — gains between -10 and 

-20 rad/ rad. The frequency at which 
this occurs may be larger or smaller 
than the basic Dutch roll frequency. 
Figure 14 shows the root locus for a 
moderate q flight condition 

(9911 N/m 2 (207 lb/ft 2 )), where 

increases the Dutch roll frequency. For 
the 0. 75L6. 1 (0. 75L20) flight condition 
(q = 18,338 N/m 2 (383 lb/ft 2 )) the Dutch 
roll goes unstable at a frequency slightly 
lower than the basic Dutch roll frequency. 
Depending on the relative locations of 
6 a 

the spiral mode pole and its associated zero on the real axis, the — 
ence spiral stability with varying degrees of effectiveness. 

6, 



(Gain values in rad/rad.) 


loop will influ - 


' a 


Figure 15 shows a typical root -locus diagram for a positive — . 


is stabilized 






The Dutch roll 
increased), and 



the spiral mode is destabilized. 
The exact shape of the Dutch roll 
locus again is subject to flight 
condition, and for some higher 
dynamic pressures the Dutch roll 
frequency increases for the posi- 
# 6a 

tive — gam. 


range of gain. 


For a realistic 
l< 10 rad/rad 


|6a 

1(3 


(Gain values in rad/rad.) 


(table, pg. 9), the Dutch roll dy- 
namics are definitely affected, 

6a 

but the — loop is not nearly as 

r 

effective in changing or 

as the ~ (figs. 8 and 10) and 
6 r P 

— (figs. 3 and 4) loops, respec- 
tively. 
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The 


Feedback Loop 


5a 

r 


u a 

The — loop commands aileron deflection proportional to yaw rate. It is not one 

of the primary control loops that would be used alone but could find applications in 
multiloop combinations or in special cases. By augmenting N r and L r , this loop 
would be expected to alter Dutch roll damping. 


jw 



(Gain values in rad/rad/sec.) 


jw 



(Gain values in rad/rad/sec.) 


A representative root -locus dia- 
gram is shown in figure 16. For the 
JetStar, the complex zeros are gen- 
erally in the right half of the plane 
for all but high dynamic pressure 
flight conditions. Relatively small 
values of feedback gain will cause 
Dutch roll instability. The spiral 
mode is stabilized as the effective 
L r is decreased. 

Reversing the sign of the feed- 
back results in the root locus of 
figure 17. The spiral mode is di- 
vergent for very small gain values, 
usually less than 0.5 rad/rad/sec 
for most flight conditions. The 
spiral instability occurs for much 
smaller gain values than does Dutch 
roll instability, thus making use of 
this loop alone to alter Dutch roll 
effects impractical. 


Jo 

The — Feedback Loop 
<P 


The — loop, commanding 


aileron proportional to bank angle, 
is somewhat distinctive in that it 
creates two aerodynamic stability 
derivatives that are nonexistent 
for the basic aircraft: 


(^) L 5 a ^’ creatin & 

(y) N 6 a ^’ creating 


L<p<P 
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The root locus for the loop is shown in figure 18. The pole-zero 

t 6 a 

configuration of figure 18 is identical to that of the - — root locus (fig. 11) except that ■ 

P 

there is no zero at the origin. The spiral mode is stabilized and the roll mode de- 
stabilized. Dutch roll dynamics are af- 
fected very little. This loop can be used 
Dutch roll V 2 to create a low-frequency mode commonly 

referred to as the lateral phugoid. In- 

5 a (^o creasing the feedback gain causes the roll 

= -1-5 — and spiral modes to merge and form an 
" i 1 oscillatory pair. The gain at breakaway 

varies with flight condition and is 
t -0. 3 rad/rad for the 0.23L0 condition. 

‘ vi spiral , „ This loop is useful in itself in stabilizing 

2 1 y m i the spiral mode or in conjunction with 

/ other loops in stabilizing the spiral mode 

— -0. 3— 1 f without interfering with the Dutch roll 

^ dynamics. 


= -1.5- 


Spiral 


5a 

— - -0. 3- 

<P 


<§ a Relating the spiral pole location to the 

Figure 18.- Root locus for the - — loop at 0.23L0. feedback gain is not easily done because 

(Gain values in rad/rad.) the spiral root cannot be simply expressed 

in terms of stability derivatives. When 
the characteristic equation of the system is written as 


As 4 + Bs 3 + Cs 2 + Ds + E = 0 


it is possible to approximate the spiral root in terms of the coefficients of equation (15). 
If the assumption is made that a single-order root exists that is much smaller than any 
of the other roots , a good approximation to this root is 


By writing the complete characteristic equation from the set of lateral -directional 
equations in appendix A, including the auxiliary feedback equation, in the form of 
E 

equation (15), the ^ term can be formed as 


_1_ 

D t s 




l A-VVvtV(* 


L )[ L 6 a (N r + 


The derivation of equation (17) is presented in appendix E. For the aircraft without any 
E 

feedback loops, ^ is usually very small. The balance between L^N r and L r N^ 
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determines the sign of the fraction and, hence, determines if the spiral root is con- 
vergent or divergent. Investigation of relative magnitudes of terms in the numerator 

6 a 

and denominator indicates that a very small variation of — is required to control 
the spiral mode. ^ 

/ 6 a \ 

In appendix E , the N™ terms, I— )N5 a , were much smaller than the L, w terms, 

/6 a \ ' ™ ' 60 ^ 

I— lL5 a , and so were neglected. Thus, the ~ loop acts as though it creates an L ^ 

alone. However, it is important to note that this does not imply that an equivalent 

will not affect the spiral mode. If can be created where Ly is insignificant, such 

6 r 

as with a — loop or with wing -tip drag surfaces, spiral stability will be greatly 
affected, as shown in reference 4. 

6 a 

The — loop induces the derivative with the conventional aileron control 

surface, through N$ a , which is a very small number for the JetStar. This feedback 

loop can also be used to destabilize the spiral mode by reversing the sign of the feed- 
back gain. The root locus for this case is shown in figure 19. The roll mode is now 

6 a 

stabilized, and the value of -y 
for which the spiral mode is 
divergent is small, usually less 
than 0.05 rad/rad. 

The expression relating 
spiral root location will give sat- 
isfactory results as long as — 

is not very large. It is also 
interesting to note that all of the 
aerodynamic stability derivatives 
augmented by previously dis- 
cussed loops appear in equa- 
tion (17). It is not surprising, 
then, that all of the feedback 
loops considered so far should 
have noticeable effects on the 
spiral mode. 




Dutch roll 

J2 

3 


- 
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Figure 19.— Root locus for the 


9 


loop at 0.23L0. 


(Cain values in rad/rad.) 


MULTILOOP RESPONSE FEEDBACK OPERATION 


The response feedback loops have thus far been considered as single-loop systems. 
In most practical applications, however, it is necessary for two or more loops to be 
used simultaneously to achieve the desired aircraft response characteristics. If, for 
example, it is required that both Dutch roll frequency and damping ratio be controlled, 
it is necessary to use more than one feedback loop, since no one loop can exercise 


It was noted that the 


independent control over both of these parameters, 
affects the Dutch roll frequency and the loop has its major effect on the Dutch roll 


ur , 

T loop 
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damping ratio, hence the two loops can be used together to adjust the system frequency 
and damping ratio. The block diagram of this system is shown in figure 20. By re- 
arranging the diagram, it is possible 

to write a (s) transfer function 

% 



for the entire closed-loop system as 


— /e\ = - 




Figure 20.- Two-loop block diagram representation. 


The characteristic equation, which is 
the denominator of equation (18) set 
equal to zero, is a function of both 
feedback gains , hence two variables 
exist which influence the movement of 
poles on the s -plane. 


The root -locus technique can be used to analyze certain multiloop systems as 
shown in reference 8. The result is a root contour, as shown in figure 21. The 


dashed lines are constant 


gains; the solid lines represent constant values of — . 

r 

iu> 



Dutch 

roll 


/ Roll 


Spiral 


Lines of constant 


■Lines of constant 


°r r 

Figure 21.— Root contour for the — ^— and ~7T~ loops at 0.23L0. 


Sr , S r 


(Gain values for — and 
respectively,) 


in rad/rad/sec and rad/ rad, 
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By using such a plot, it is possible to achieve a desired c and £w,. The roll 

5 r 6r 

.and spiral modes are not shown in the figure, but each combination of — and — 

will determine a roll and a spiral time constant. If a somewhat independent means of 

6 a 

varying T r is required, it will be necessary to add another loop, — , for example. 

P 

Spiral stability may also be important; hence, a fourth loop would be necessary. Thus, 
there are many combinations of loops which could be used, and these configurations 
are dictated by the requirements on the closed-loop response. 


The method for analyzing a multiloop feedback system involving a single control 
input, 6 a or S r , can be summarized as follows: 

r r 

1. Close one loop and find the root locus corresponding to various values of 
feedback gain. 


2. For one value of the feedback gain of loop number one, locate the roots 
(poles) and add the zeros of loop number two, disregarding the original zeros. 

3. Find the root locus corresponding to the pole-zero configuration of 
step 2 , as a function of the second feedback gain. 


Note that step 3 may be repeated for every value of the first feedback gain desired, in 
each case using the pole locations corresponding to that gain. The resulting diagram 
is a root contour, showing graphically the effect of two feedback loops. The process 
can be expanded to more loops only if the loops have the same control -surface input. 

6 r 6a 

This method is not valid for a — and a — loop combination, for example, inasmucl 

as the implicit assumption of this technique is that closing the first loop did not alter 
the zeros of the second loop in figure 20. This assumption follows from the general 
statements made earlier concerning the effects of feedback loops on transfer-function 
numerators. 


For multiple input-output systems , the characteristic equation must be written 
for the total system and factored to yield the roots which determine the transient 
response of the aircraft. This must be done for each combination of feedback gains. 
Reference 9 presents a more detailed analysis of general multiloop systems. 

In some cases, it is possible to develop an expression relating two feedback gains 

and a dynamic mode variable. Equation (19) applies when the 
used together as 

6 r _ 2 + N r + ^ ^ 

p ~ 


— and — loops are 

0 P 


(19) 
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Equation (19) was obtained by examining equation (6) and noting that 
Y p and indicates the required — gain for a desired 2£^o>^. The 
which results from the loop used alone. 


(f) 


CO,, 


Ys r augments 
now is that 


In most cases it is not possible to write simple expressions relating two or more 
feedback gains, and it is necessary to resort to fectoring the total system character- 
istic equation, given the feedback gain values, to find the closed-loop dynamics. 


ANALYSIS OF TWO-LOOP OPERATION FROM 
STABILITY-DERIVATIVE CONSIDERATIONS ALONE 


In some cases , the use of feedback loops to alter the dynamics of an aircraft is 
not conveniently analyzed by root -locus methods. Although quantities such as frequency, 
damping, and time constants are easily picked off s -plane root -locus plots, other 

parameters often used in handling -qualities studies are not, such as or |^-(s)| . 

Often, these parameters are of primary concern in a response feedback loop appli- 
cation. Analysis of such applications is possible by working with the equivalent aero- 
dynamic stability derivatives, that is, the effective stability derivatives, resulting 
from feedback loop augmentation. 


In many handling -qualities studies, it is desirable to vary the ratio £(s) , 

Ip 

Reference 2 contains the following simplified expression for this ratio when tip is 
small: 


£ 


(s) 


L£ 

N/3 


effective 


( 20 ) 


Any feedback loop which augments L^ or N^ will affect the (s)j^ ratio of equa- 
tion (20). Rewriting the equation to include the effects of a 


6 r , \ ^ 6 r 

— loop (where 

carries its own sign) introduces new terms as shown in the following expression: 


£ 


(s) = 


L0+ 

8$ 

L5 r 

N^ + l 

U ) 

K 


( 21 ) 


Changes in the 


^(s) ratio are necessarily accompanied by changes in the Dutch roll 

Ip 1^ 5 a 

natural frequency because of the N^ augmentation. The use of alone is likewise 
undesirable because of its influence on the effective N^ through the term ^pyNg a , 
as shown by the following equation: 
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By using an expression for 


£ 


(s) 


which includes the effect of both the -jr- and the a 


, ... ll P \(p I P p 

loops, it is possible to obtain a desired kg ( s )| ratio and yet maintain the effective 
constant, assuming that is an accurate approximation to o^ 2 , that is, the o; t L^ 
is small due to a small a^. This expression is 


P 


IH 


J /3 + 






N/3 + 


(tMt) 


N S fl 


(23) 


If the total effective is to remain constant at the open-loop value, the relation 


6r 

P 


a\ N ^a 

V I N6 r 


(24) 


must be satisfied. Substituting equation (24) into equation (23) results in the expression 



(25) 


The desired 


£ 


0 (s) ratio can be obtained by solving for the — gain in equation (25), 

P 5 r p 

then solving for the corresponding gain from equation (24). The absolute -value 
sign indicates that there will be two sets of gains which will satisfy equation (25). The 
selection of the sign of the ^ gain to be used must be made in the light of other con- 
siderations, such as the fact that is an important handling -qualities parameter in 

itself, and for conventional aircraft in the transport category has a minus sign 
associated with it, which indicates positive dihedral. This means that the sign of the 

— gain must be chosen that will maintain the total effective as a negative quantity. 


Other dynamic side effects of the two -loop combination used here, such as the 
influence on the spiral mode, may dictate the use of other loops to maintain a flyable 

vehicle. The effect of these other loops on the ms) ratio, in turn, may have to be 

\P \ip 

taken into account. Iterative processes will usually be necessary to determine the final 



combination of gains, and these processes are best handled by a digital computer. For 
some cases when equation (20) is valid, equation (25) may be adequate to calculate the 
required gains. In the preceding example, it was assumed that it was necessary to 
hold N p constant at the open -loop value, but equation (24) could just have easily been • 

written to attain some other desired effective N^. In practice, there will be many 

cases in which equation (20) will not be adequate. This must be determined first. The ‘ 

simplified approximation to |^(s)| was used to illustrate the equivalent stability - 

I P |i p 

derivative approach in a two-loop situation. 


EFFECT OF ACTUATOR DYNAMICS ON RFS PERFORMANCE 


Analysis of the response feedback loops has been carried out on an idealized 
system where neither sensors nor actuators display any phase lag. The scope of 
this paper does not warrant a detailed examination of all the associated effects of 
sensor and servo dynamics; however, it is worthwhile at this point to consider how 
these effects can be included in the analysis of a typical feedback loop. 

The problem of phase lags in a response feedback system is a serious one. 
Generally, sensors, such as gyros, accelerometers, and vanes, can be expected to 
have less phase lag than the actuators at operating frequencies. The effect of actuator 
dynamics on response -feedback -system performance is, then, of primary interest. 

A suitable mathematical model must be chosen to represent the actuator. One 
such model is shown in figure 22. The actuator can be shown in the aircraft response 

feedback system in the forward loop 



Figure 23.— Inclusion of actuator in RFS loop. 
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The order of the characteristic equation is increased from a fourth- to a sixth-order 
equation. The phase lag of the actuator itself can be calculated at any driving frequency 
as 



(27) 


It is apparent from equation (27) that in order to minimize phase lag co n should be as 
large as possible for a given £ n . 


The addition of an actuator in the forward loop does not greatly increase the labor 
required for analysis by the root -locus method. The poles (and zeros, if any) of the 
actuator transfer function are placed on the s -plane with the poles and zeros of the 

§—( s ) transfer function. The location of the actuator poles is a function of £ n and co n . 

It is usually desirable to have approximately 70 percent critical damping for the 

actuators. The effect of such a 70 -percent damped actuator on the loop and its 
influence on the Dutch roll are shown in figure 24. ^ 



Figure 24.— Effect of actuator on Dutch roll dynamics with 
8 r 

the — loop closed at 0.23L0. 

P 


The locus for the system when 
actuator dynamics are neglected 
was discussed in the section on the 



has destabilized the Dutch roll by 
decreasing the damping ratio. For 
an actuator with a natural frequency 
of 25 rad/sec, the Dutch roll can be 
made to go divergent if the gain is 
high enough, in this case, above 
20 rad/rad. Actuators with natural 
frequencies less than 25 rad/sec 
would cause instability at lower 
values of feedback gain. 

The roll and spiral modes will 
not be affected noticeably unless 
they form a lateral phugoid and 
break away from the real axis. 

The breakaway point itself will be 
dependent somewhat on actuator 
pole location. The representation 
of the actuator may be different 
from the ideal second-order type 
assumed here. Position rate or 
pressure rate feedbacks around the 
actuator, along with other compen- 
sation, could result in a higher 
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order representation with possible numerator roots. In this case, the roll and spiral 
modes would be influenced. It is also important that the actuator poles be examined 
as the feedback gain is increased. While it is possible that these poles will move into 
the right -half plane, resulting in an unstable system, it is also likely that they may 
move toward the real axis, indicating that the damping ratio of the actuator is in- 
creasing. Equation (27) shows that the phase lag will increase with increasing £ n , 
and it is undesirable to add any more phase lag to the system. 

The importance of including actuator dynamics in practical applications is evident 
from figure 24, and especially when high gain feedback loops are used with low fre- 
quency actuators. 


CONCLUSIONS 


On the basis of an analysis of the response feedback system for the lateral- 
directional modes of a variable -stability transport airplane, the following general con- 
clusions can be drawn: 

1. Each feedback loop influences all of the aircraft dynamic modes but, in most 
cases, displays one dominant effect. 

2. For values of feedback gain that are considered to be reasonable for this air- 
plane, most classical approximations of loop action were valid. 

3. The analytically derived gain expressions have conditions of validity associated 
with them, and it is important to ascertain that these conditions are satisfied. 

4. The equivalent stability -derivative method of analysis is valuable when root- 
locus diagrams are not convenient to show loop effects. 

5. Root -locus analysis techniques can be expanded to miltiloop configurations, but 
it must be realized that the methods outlined herein are applicable only to loop com- 
binations that do not alter the numerators of transfer functions in the closed loop. 

6. Actuator dynamics can be readily included in a root -locus analysis. In more 
detailed studies of response feedback systems it will usually be necessary to include 
actuator effects. 

7. Practical limitations will determine the final success of a response feedback 
system. These limitations must be included in the analysis as demanded by the 
particular application. 


Flight Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif. , December 14, 1966 
720-04-00-01-24 
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EQUATIONS OF MOTION USED IN LATE ML -DIRECTIONAL RFS STUDY 


The side -force and moment equations (ref. 11) are written in a set of orthogonal 
body -axis coordinates. The three equations written in conventional nonlinear form are 
as follows: 

Side force - 


• 1 9 

m(v + ru - pw) = mg sin (pcosQ + ^PV S(-C D sin/3+ Cycos/3) 


Rolling moment - 

p j xx - rI xz + ( J zz - !yy) c i r - p^xz = |pv 2 sbc z 


Yawing moment - 


rlzz " P J XZ + ( J YY - x xx) PPl + *l rI XZ “ 


|pV 2 SbC n 


These three equations are modified by assuming that, in still air, 



P 




To further reduce complexity with little loss in accuracy, the following approximations 
are made: 

1. Inertia coupling terms in the moment equations are dropped 



3. -Cj) sin /3 + Cy cos j3 is approximated by Cy 
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4. V is used for u in the side -force equation 
The side -force equation then becomes 


Dividing by mV, 


• 1 o 

m(v + rV - pw) = mg sin <p cos 0 + — pV^SCy 


v w g* oVS 

y + r - Py = y sin 9 cos 9 + 2^T C Y 


Incorporating the approximations , 


b = ' 2 ^’ C Y + y sin c °s © + 


pa - r 


The rolling -moment equation is 


-i— /iT„„ j. ir>v2< 


P = j^XZ + g-pV^SbCj) 


The yawing -moment equation is 


* - f^XZ + p v2 SbC n ) 


The three equations are then linearized about a trim flight condition (V = Vy) where 
the initial p, r, <p, i p, p, 6 r , and 6 a are zero. The resultant linearized 
perturbation equations are (assuming Cy^ = Cy^ = = 0) 

Side force - 

, VyS / \ VrpS / u \ „ 

P=|J S (°V + C *6 T 6 rj+ P"25 (pvj) C *r r + V* P“T ' r 


Rolling moment - 


i> = 5 ' + S PVT2Sb ("V + V r + Ci «a 6a ) + 2 k r + % P )] 
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Yawing moment - 


r = ^ P+ ^[^ pV T 2sb ( C V l+ °»6 r 6r + C »a a 8a ) + l ,,V T 2sb ( 2 vi)(Cnr r+ <V>)] 


As a further modification to conserve analog equipment necessary for programing, the 
Euler angular rate relationship is solved for p as 


p = <p - ip sin© 

Substituting the Euler expression for ip, 

P = <P ~ (q si n <P + r cos <p) tan 0 

Linearizing this expression and assuming tan 0 T = tan a T = o; T , 


p = <p - ra-p 


Dif fe rentiating , 


p -ip - ra T 


Substituting p and p into the three equations and rearranging, using the dimensional 
forms of the aerodynamic stability derivatives, yields the final form of the lateral - 
directional equations of motion. Note that the equations are now in terms of an earth 
reference, <p and ip, rather than a body axis, p and p, as follows: 

Side force - 


(1 - Y r + a T 2 )r + 4 - - a T <p - ^cp = Y 6r 6 r 

Rolling moment - 

-(a T + r + (a T L p - L r )r - L^+ <p - L p y = LQ a 6 a + L 6r .6 r 

Yawing moment - 

(l + Of T 1 ~) r + (a T N p - N r )r - '<P - N p ^> = N 6j 5 a + 


r 
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By Laplace transforming the equations and using the variable s, the equations can be 
written as 

Side force - 

(1 - Y r + o; T 2 ) r + (s - Y^/3 - (Q! T s + ^-) <p = Yg^ 

Rolling moment - 

[-(«T + s + < a T L P " L r)j r - L^/3 + (s 2 - L p s) <p = Lg a 6 & + Lg r 6 r 

Yawing moment - 

K 1 + “T m) s + (“ T N P - N r )] r - s2 + Np s) V = N 6a a a + Nj r 6 r 
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APPLICATION OF ROOT -LOCUS METHODS TO 
ANALYSIS OF RESPONSE FEEDBACK LOOPS 


The root -locus method of analysis (refs. 7 and 8) directly relates feedback gain 
values and characteristic equation roots. It is these roots that determine the form of 
the response of the dynamic system. Generally, the characteristic equation for the 
lateral -directional mode is a fourth-order equation which, for the transport class of 
aircraft, usually factors into the following form (ref. 5): 

(s 2 + 2fyotys + oty 2 )^s + ^-)|s + = ° 

Each parenthetical term is referred to as a dynamic mode. The second-order 
term is the Dutch roll mode, and the two first-order terms are the roll and spiral 
modes, respectively. Adequate descriptions of the physical characteristics of these 
modes are found in most texts on aircraft dynamics. For the analyses of this paper it 
is sufficient to note that each parenthetical term in the equation in the Laplace variable 
(frequency domain) has an equivalent form in the time domain, as explained in refer- 
ences 7 and 8. The Dutch roll is an exponentially varying sinusoidal oscillation of the 
general form 


where A is an arbitrary phase angle. The roll and spiral modes are, respectively, 

_ t _i_ 

exponential functions of the form e and e T s . The root -locus technique provides 
a convenient method of plotting the roots of the characteristic equation as a function of 
the particular feedback gain used. The resultant curve, presented on the s -plane, 
yields information about the dynamic stability of each of the modes. Figure 25 shows 
the three modes and indicates the parameters characteristic of the dynamics of each 
mode. 

The basic aircraft characteristic equation roots are designated by X's and are 
referred to as poles. The Dutch roll mode appears as two poles because the roots 
of the second-order term in the characteristic equation are a complex pair. The 
magnitude of the distance between the origin and one Dutch roll pole is the undamped 
natural frequency in radians/second. The component of this distance along the 
imaginary axis is the damped frequency that would be measured on a time history. 

The damping ratio of the Dutch roll is the cosine of the angle between the negative real 
axis and a line drawn from the origin to the pole. The roll and spiral mode time con- 
stants, in seconds, correspond to the inverse of the distance between the origin and 
the pole location. If a root lies to the left of the jco axis, the resultant time response 
is stable; whereas, roots falling to the right of the jco axis are unstable. Roots 
falling on the jw axis are termed neutrally stable. 
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Graphical techniques used in plotting root-locus diagrams utilize the transfer 
function numerator roots (zeros) as well as the poles. These zero locations are 
designated by O’s and are shown herein on root -locus diagrams. The graphical 
techniques are described in references 7 and 8. 


Imaginary axis 



X 

Dutch roll 

Figure 25.— Lateral-directional modes as displayed on the 
complex or s-plane. 
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TWO-DEGREE-OF-FREEDOM APPROXIMATION TO 

5 r 

THE DUTCH ROLL WITH THE — LOOP 

r 

By choosing only the dominant terms in the equation of motion, the side -force and 
yawing -moment equations become 

(s - Yp)0+ r = Y 6 r 6 r 

and 

-N 0 / 3 + (s - N r )r = N 6 r 6 r 

allowing only (3 and r motions. Now letting 6 r = 6rp + 

(s -Y^)/3+r = Y 6 r 6 rp + ^Y 6 r r 

-N^/3 + (s - N r )r = N 5 r 6 rp + (y)N 6 r r 
or 

(s - Y0P* [l - ^f)v 6 r ]r - Y 6 r 6 rp 

-Xffi * [* - Nr - N6 r ] r - N 6r 5r p 

the characteristic equation is 

s 2 - [n r + Yp + (^)n 6i ]s + Np + N r Y ? + (^)(NS r Yp - N„Y 6r ) = 0 

In similar fashion, writing 5 r = 6 r the two-degree-of-freedom approximation 

to the Dutch roll can be used to obtain the characteristic equation 

s 2 - [Nr + Y „ + (j) Y 6 r ] * + V * N r Y^ + J K * N r Y 6r ) = 0 
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ANALYTICAL EXPRESSION RELATING 
2 fyoty OF THE DUTCH ROLL 


«r 

P 


TO 


The expanded characteristic equation general literal form is easily shown (ref. 5) 
to be 


s^ + 





+ 



s 


3 


L 2 + + _J_\ 

\ V T r T s r r r s J 


+ 




+ 2 i ^ 

T r T s/ 


s 



Comparing the coefficients of this equation with 


s 4 + | s 3 + S s 2 + P 


s + 


I- 


which is the general form of a fourth-order algebraic equation where the coefficient of 
the highest order term s^ has been made equal to unity, and equating the s 3 terms 
of the two equations results in 



ip^ip + 


1_ 

T r 


+ 


1 


T s 


For most transport aircraft, — can be neglected with respect to the other two terms. 

T s 

Also, to a very good approximation 


and so 



A 2 h LOl P L 1 
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After writing the feedback loop equation as 

and substituting into the lateral-directional equations of motion, the characteristic 
equation can be found. The coefficient of s 4 term (A) is 


I 1 -(xHK i+ (“ t 



o 

and the coefficient of s term (B) is 


l 1 '(r) Y6 J[ ( “T N p - Kr > + iff <“T L p - L r>] 
-( 1+ “ T lfi)[ L P + Y /3 + (fj f T L «r - L P Y 6 r )j 

4 


Ixz 

Ixx 


j| :, » - V «H 




The s 4 and s^ coefficients can be simplified by neglecting small triple -product 


XZ 


small numbers themselves, appear in 


ixz 

terms; inertia terms, since -7 and — J 0111UXA UUiIlkFVi.0 U1W111UV1VV/D) X XXX 

ixx -■‘Z z 

products usually much smaller than other terms in the sum; and terms involving a^, 

since the trim angle of attack is assumed to be small. Using the remaining significant 
terms, the first two terms of the characteristic equation are 


and 



But, as shown previously, 


B 

A 


~ kp 
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and, substituting the expressions derived for A and B, 


B 

A 





Equating the two expressions for 


B 

A 


and solving for 



9 


6 r + N r + Y /3 

/3 ^6 r + Y 6 r |^r + 
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* 


« 


APPROXIMATION TO THE SPIRAL MODE TIME CONSTANT 


6 


a 


WHEN -=■ IS USED AS A FEEDBACK LOOP 


Reference 5 states that if the spiral mode constant r s is very large, 


1 


T 


S 


ca 


E 

D 


where D and E are the coefficients of the s 1 and s° terms, respectively, of the 
lateral-directional characteristic equation. If the auxiliary feedback equation is 
written as 


«a = S + (^V 

D, the coefficient of the s 1 term, is 


D — -(oirjiNp - Nj.) 


+ (Q!.pLp - L r ) 


[“T L /i - L p Y /3 + (^ 

jK] 

[«T N ,! - N p Y ,i + (- 

iK] 


-<1-Y r + a T 2 ><W - L ^ N p) 


The E term is 


E = -(« T N p - N r )[^ L 0 -(^^j 
+(Of T Lp - L r ) |y^ ~(^ N fi a Y ^J 
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(1- Y r +Q! T 2 )B L 6 a ^-(f)^N 6 


Neglecting small triple -product terms and assuming that inertia terms are small, 

ixz 

since - — « 1 for the JetStar, 

!ZZ 

D - L^Np - L pN( j - ^ L(J + (|)[L 6a (N r + Yp) + N 6a (a T Y„ - L^] 


For the JetStar, 


L6 a (N r + Yp> »N 6a (a T Y S - L r ) 


D = L /3 N p - L p N /3 * Vt: L 0 + 


(^)k^ + Y o)] 


Similarly for the E term, neglecting small triple-product terms and inertia terms, 
E = f T (¥r - E rN 3 ) -^)[L6a^rY^ N„) - N 6a (L r Y^ + L„)] 


For the JetStar, 


N/3 »N r Y^ 


» L r Y^ 


E = #:( L /3 N r - Wfi) -^)(L 6a N^ - ^N 6a ) 


Furthermore, for most cases, 


L 6 a N /3 » l /3 N 5 £ 


and so the fraction — can be written as 


E M J_ 

D r„ 


^(y<r - LrNp) 

Vp - LpNp - ^ +(f)[L« a (N r+ Y„) 
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